Dealing with delicate issues in waveforms calculations 
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We revisit the calculation of gravitational radiation through the use of Weyl scalars. We point out 
several possible problems arising from gauge and tetrad ambiguities and ways to address them. Our 
analysis indicates how, relatively simple corrections can be introduced to remove these ambiguities. 
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I. INTRODUCTION 



The definition of gravitational radiation involves a 
fair amount of work so as to distinguish variations in 
the metric that are purely coordinate effects from true 
signals propagating on a background which itself must 
be well established. At finite (but sufficiently far) dis- 
tances from an isolated source several approaches ex- 
ist to do so provided a suitable background can be 
identified [3, 0, E II H • For sufficiently generic cases the 
unambiguous identification of the gravitational radiation 
can only be done at future null infinity (I + ). This is pos- 
sible since when studying an isolated system, future null 
infinity is the asymptotic region of an asymptotically flat 
spacetime. There the metric is exactly flat and distur- 
bances around it can be associated with outgoing gravi- 
tational waves subject to singling out a preferred frame. 
This issue arises as there is no unique flat asymptotic 
metric. In fact there are as many of them as supertrans- 
lations are in the BMS group; the asymptotic symmetry 
group. Nevertheless there is a well defined notion of grav- 
itational radiation which exploits the algebraic notion of 
asymptotic flatness [6] which in turn implies the existence 
of the gravitational radiation fields ^ and ^3. The ex- 
istence of these fields can also be deduced from stronger 
conditions, like peeling]?], HJ. Then, after a careful choice 
of coordinates, radiation can be unambiguously defined. 

Unfortunately, numerical applications dealing with 
black hole spacetimes can not yet reach I + un- 
less Cauchy-characteristic extraction is employed^, [1(3] 
for the extraction procedure or Cauchy-characteristic 
matching[Tl| or the conformal formulations [l2| are fully 
realized to model the whole spacetime. Without adopt- 
ing any of these options a commonly employed approach 
for obtaining gravitational waves relies on applying the 
infrastructure developed for X + at finite distances. Here, 
several compromises are made as natural quantities de- 
fined at T + need to be translated somehow to finite dis- 
tances where they need not be well defined. 

The purpose of this note is to point out difficulties that 
likely arise in this approach and a route to address them. 



II. REVIEW OF THE APPROACH 

The calculation of gravitational radiation at I + is 
based on the foundational works of Bondi[l3[ and 



Sachs LJ| as an approximation expansion at infinite dis- 
tances and by the work of Penrose [lj| as a geometrical 
construct at the boundary of the physical spacetime. Un- 
der reasonable assumptions one can show that the space- 
time is asymptotically flat and time dependent perturba- 
tions close to the boundary are either due to gravitational 
waves or coordinate effects. In order to disentangle coor- 
dinate with physical effects, care must be taken in dealing 
with the asymptotic symmetry group and fixing a conve- 
nient frame with respect to which inertial observers can 
define radiation. 

Our goal here is not to revisit this approach (of which 
descriptions have been presented in e.g. [8l. [iH \vK [iH [iH 
|20J]). but simply to state the crucial ingredients which are 
missing when the strategy is employed at finite distances. 
In the next sections we will carry out an analysis of what 
corrections are required when key assumptions described 
below are missing. 

These ingredients are, 

• Peeling is assumed, this means that suitable Weyl 
curvature components behave asymptotically in a 
well defined manner. 

• Outgoing null hypersurfaces, parameterized by u 
intersect I + defining a sequence of S 2 surfaces. 

• Each of these surfaces is conformal to a unit sphere 
metric; off this surface (into the spacetime) the de- 
parture from it is of lower order. Namely the an- 
gular metric in a neighborhood of X + can be ex- 
pressed, as g AB = r 2 h AB = r 2 (q AB /V 2 + c AB /r + 
0(r~ 2 )); with V a conformal factor, q AB a unit 
sphere metric, and r a suitably defined radial dis- 
tance. 

• A null-tetrad {£ a , n a , m a , rh a } satisfying £ a n a = 
~m a rh a = 1 (with all other products being 0). 
The spacetime metric can be expressed as g a b — 
2^(a n b) — 2TO( a m{,), and suitable projections of the 
Weyl tensor can be defined. 
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Using standard *2l],[22| conventions for the Riemann 
tensor and spinor dyad, particularly useful scalars 
obtained from them are, 



*4 = C abcd n a m b n c m d 
*3 = C abcd e a n b fh c n d 
* 2 = C abcd e a m b m c n d 



(1) 
(2) 
(3) 
(4) 



• A suitable (Bondi type) expansion in terms of 1/r 
coupled to coordinates chosen such that V = 1, 
g® r = 1 and g° A = (x A labeling angular coordi- 
nates at u = const, r — > oo) gives rise to several 
important relations @, Q . In particular, 



-3cr 
9*° + a°¥ 4 



(5) 
(6) 
(7) 



where the supra-index "°" indicates leading order 
in an expansion in the radial coordinate and 3 is 
the edth operator[2lj of the unit sphere. 

With this structure at hand, the following results are 
obtained. Given any section S at future null infinity the 
Bondi momentum is given by 



P a = -— [ i a (V° 2 + a°a°)dS 2 
4?r J s 



(8) 



where l a — (1, sin(0) cos(^), sin(0) sm(^), cos(0)) when 
expressed in standard angular coordinates (9, <f>) . The 
so called Bondi mass M is the timelike component of 
this vector; namely 



M 
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(9) 



It is interesting to note that the Bondi momen- 
tum can also be expressed in the terms of the Psi 
supermomentum[23| ^ = + <j°a° + 3 2 ct° by 



P a = / /"*r/,S' 

47T 



(10) 



since it has a couple of useful properties; namely it is 
real, '& = and its time derivative is simply 



(11) 



Therefore, the time variation of the Bondi momentum, 
in terms of time Bondi coordinate, is just given by 



pa _ _ 



1 

4tt 



i a &°a°dS 2 . 



(12) 



Notice that due to relation |5]) one might choose to 
replace any appearance of &° by the time integral of ^9. 
This is commonly done in numerical simulations [2J, |25| . 
S H3, S S IS, EH, [S3, S3 ; however one should notice 
that is only valid in terms of a Bondi tetrad and 
coordinate system. 



III. CONNECTION WITH FINITE DISTANCES 

The above expressions provide a well defined infras- 
tructure with which gravitational radiation can be de- 
fined. However, to obtain them one has introduced key 
ingredients, otherwise the expressions would result sig- 
nificantly more involved. In particular without choosing 
the conformal frame and/or radial distance such that the 
leading part of the angular metric is indeed that of the 
unit sphere, several changes arise even when all other 
properties are satisfied. For instance, if the angular met- 
ric were given by $ab = r 2 (qAB /V(u,9,<p) 2 + 0(r -1 )) 
then, instead of §5§ one would have 



1 



2 • 



(13) 



Therefore, if we call the inertial (Bondi) radiation 
field one would have the relation 



1 ( ^ 3 • 




where the factor V~ 3 results from one of the corrections 
which we discuss in this work. Notice that this is not 
a perverse circumstance, rather it can be generically ex- 
pected as surfaces at u = const, r = const define a 2- 
sphere whose metric is always conformally related to the 
unit sphere. 

Certainly, in principle coordinate conditions can be 
adopted so that at the extraction sphere the simplifying 
conditions hold, these conditions need not be those that 
are convenient for the numerical implementations. In- 
deed, coordinate conditions are exploited in a non-trivial 
manner to aid in the numerical simulation [3^. l35l |36| . 
Therefore, one is then left with having to consider how 
to correctly identify suitable coordinates and/or a frame 
to extract the desired quantity. Here, one of two op- 
tions can be adopted. One approach is motivated from 
a perturbative point of view and relies in extracting a 
suitable Kinnersley tetra d 13711 from the numerically gen- 
erated spacetimefH, [39|, l4dl |. The so called Kinnersley 
tetrads are null tetrads that were constructed for the 
study of type D vacuum spacetimes; but one can in gen- 
eral choose them so that to leading order they conform 
to a Bondi tetrad. The basic idea is to use a null tetrad 
adapted to the asymptotic structure so that in the limit 
coincides with a Bondi tetrad. Armed with this tetrad 
the computation of the radiative properties of the space- 
time is then well defined. The delicate point however 
stems from the fact that this approach requires a suit- 
able numerical identification of a Bondi frame where er- 
rors can arise. The influence of these errors in a related 
context has been recently discussed infill]. 
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A different approach is to consider the asymptotic 
structure of asymptotically flat spacetimes and define 
the analog quantities at finite distances. This approach, 
which we adopt here, does not rely on obtaining back- 
ground quantities correctly, rather one needs only iden- 
tify a suitable Bondi frame which will essentially be un- 
ambiguous as long as the extraction worldtube is suffi- 
ciently far. We will not concern ourselves with how far 
must the extraction worldtube be, rather we assume one 
is able to place it sufficiently far for the peeling property, 
in the extracted quantities, of the Weyl scalars to be ob- 
served. Under this assumption we analyze what conse- 
quences, and most importantly, modifications one must 
take into account when calculating physical observables 
with the coordinates naturally induced on the extraction 
worldtube in a numerical simulation. 

To this end, it will be important to calculate the in- 
duced metric on the extraction worldtube and evaluate 
different quantities which play a key role. We begin by 
following precisely the same procedure commonly em- 
ployed in numerical simulations (see for instance p^. [28l . 
[42[). First a Cartesian extraction worldtube is defined 
by x 2 + y 2 + z 2 = R 2 . Observers at this worldtube are 
defined by their trajectories on the worldtube given by 
(t, x = const., y = const., z = const.). Then, 'J 4 is calcu- 
lated on the hypersurface E t (defined by t = const.) by 
suitably adopting a tetrad frame. This tetrad is chosen in 
a straightforward manner by adopting three orthonormal 
spatial vectors {r a , 9 a ,4> a } on E t (which are the analogs 
to the radial and tangents vectors to the sphere) and 
combine them with the unit normal vector N a to the hy- 
persurface. Once ^4 is obtained, an interpolation on to 
the extraction sphere is carried out. At this point it is 
useful to recall that the extraction must be done along 
null rays, thus any comparison (or refinements of the 
extracted quantities [49]) at different worldtubes should 
contemplate suitable time offsets. This offset reflects the 
arrival time of a null signal from a given internal world- 
tube to the next. 

As mentioned, this procedure, need not give rise to a 
structure fully compatible with those mentioned in sec- 
tion I. In particular the induced angular sphere metric 
need not be sufficiently close to that of the unit sphere, 
inertial observers coordinates might 'shift' in time around 
the worldtube and their associated times might tick at 
different rates. At finite distances these issues are even 
more relevant since the waves themselves will influence 
the geometry of the extraction worldtube as they prop- 
agate across it. In the next section we discuss the cor- 
rections required and will illustrate their application in 
section V. As we will see, even in linearized problems 
these issues can play a non-trivial role. 

A short summary of what these corrections entail is (in 
addition to the standard first-step to calculate ^4) is, 

• Assuming that the extraction world tube has been 
chosen far enough [HfJ, calculate the leading order 
behavior vf^ from ^4 = ^"r -1 +0(r~ 2 ); where r is 



an appropriate asymptotic radial coordinate; such 
that r = R at the extraction world tube. 

• In a similar way, obtain the leading order behavior 
of the induced metric gr on the sphere at the ex- 
traction worldtube defined by t — const., r = R, 
from gr = —qR 2 V~ 2 + 0{r r )\ where q is the unit 
sphere metric. 

• Obtain the leading order behavior g® r from the ex- 
pansion g ur — g\ r + 0(r _1 ); which measures the 
time observers of ^4 have. 

• Obtain the leading order behavior g^ A from the 
expansion g u A = g® A + C(r _1 ); which measures 
the coordinate shift in time observers of \&4 have 
around the worldtube. 

• Bondi's radiation ^ is then obtained by relatively 
simple correction factors from the knowledge of 
{j2 r , 5°^, V}. In the particular case where g® A ~ 
the expression is 



IV. ASYMPTOTIC STRUCTURE OF 
ASYMPTOTICALLY FLAT SPACETIMES 

In this section we examine which corrections must be 
taken into account to remove the ambiguities described 
previously. To this end, we will examine how the rele- 
vant quantities transform among different coordinate sys- 
tems compatible with the asymptotic structure. Then, by 
adopting one of these systems to be a Bondi one, we will 
be able to extract the correcting terms. For the sake of 
simplicity in the derivation we adopt the spinor-calculus 
approach, though the results obtained are completely in- 
dependent of this technique. 

We begin by determining key relations between differ- 
ent frames, to do so it is convenient to have at hand some 
basic equations of the asymptotic structure. 

A. Basic variables 

The asymptotic geometry can be expressed in terms of 
a complex null tetrad (£ a , m a , m a , n a ) with the proper- 
ties: 

g ab l a n b = -g ab m a m b = l (16) 

and all other possible scalar products being zero, the met- 
ric can be expressed by 

g a b = l a nb + n a £ b ~m a m b - M a m b . (17) 

Such a null tetrad is easily related to a dyad of spinors 
(o a ,l a ). The relation of the null tetrad with a spinor 
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A L A ' 



dyad is given by i a o o , m a < 

fh a -i=> i A o A and n a l a l a . Thus, determining 

how the tetrad trasnforms one infers the transformation 
of (o a ,l a ). In our discussion, to make a direct con- 
tact with the standard literature on the subject we also 
adopt a null polar coordinate system (x° , x 1 , x 2 , x 3 ) — 
(u, r, (C + C)i j{C — C))- However, as opposed to the 
standard discussions we will not make use of the avail- 
able coordinate freedom to simplify the treatment at the 
onset. Rather, we will adopt coordinates consistent with 
those employed in the extraction procedure at finite dis- 
tances and deal with the consequences of this choice. 

In most numerical applications, one normally works 
with finite size regions; and would like to estimate the 
asymptotic fields in terms of null tetrads based on a 
choice of coordinate system. At the extraction world- 
tube, one assumes that a coordinate system (t, r, C) 
can be constructed. Let L, be the timelike surface de- 
fined by (r = R = const.). On L, a null tetrad £, n, to, to 
can be defined in the following way. First, define the null 
function it such that on T one has u = t; and I = du 
everywhere. The function u is chosen such that the fu- 
ture directed vector £ points outwards with respect to the 
2-surfaces (t = const., r = R); which are, topologically, 
two dimensional spheres. 

The complex vectors m and to are defined to be tan- 
gent to the spheres (t — const., r = const.). Furthermore, 
one can choose to to be proportional to ; and to to be 

proportional to in the asymptotic region for large r. 

Then, by requiring the standard normalization condi- 
tions given in (|16p one settles the remaining null vector 
n a . 

The coordinate system (u, r, (, £) is thus straightfor- 
wardly related to the one commonly employed in numer- 
ical efforts (t, R, 6, <f>) though we will continue our dis- 
cussion with the former as it is the one employed in the 
standard literature on the subject. The conclusion how- 
ever, will be independent of this choice. 

Keeping (( = const.) and (( = const.) on the null hy- 
persurface u = const., and increasing r one moves along 
a null direction. Since I is contained on the hypersur- 
face u — const., one deduces that £ is proportional to J^. 
Then, one can write 



\9ur or 



(18) 



The appearance of g ur in the denominator is the first cor- 
rection to be accounted for due to the structure induced 
at r not necessarily conforming to a Bondi frame. The in- 
tersection of the null hypersurface u = const, with future 
null infinity defines a two dimensional surface denoted by 
S. A natural null tetrad is then completed with, 



4 = (du) a 



1 dx l 



(19) 



(20) 



" ' dx i 



du J \dr ) \dx l 



with i = 2, 3 and components £*, U and X % are: 



? = k + o 



f 3 ( 1 



(21) 



(22) 



(23) 



with 



(24) 



where V = V(u, £, C) and the square of Pq = 
the conformal factor of the unit sphere; 



_ (i+CC) 



is 



U = rUoo + U + — + O 



1 



(25) 



where 



Uoo = U = —AV, U x 



2 - 2 , (26) 



where Ky is the Gaussian curvature, given by 

K v = 2V 99 V - 2dV dV + V 2 
= VV 2 V - V'VViV + V 2 
= V 2 V 2 \nV + V 2 , 



(27) 



of the 2-mctric 



dS 2 



V 2 P 2 



d( d(; 



(28) 



where the regular conformal metric restricted to S is pre- 
cisely g \s= —dS 2 , V 2 = V l Vi is the Laplacian operator 
of the unit sphere metric and Vj its covariant deriva- 
tive. Let us emphasize that V = 1 makes dS 2 in (|28|) the 
unit sphere metric. Finally, the other components of the 
vector n a have the asymptotic form 

A 2 = A 2 » + (-L), X 3 =X! + 0^y (29) 



B. Asymptotic gauge freedom, restricted case 

At this point, we find it convenient to restrict to a sim- 
plified case before describing the general one. We will 
assume here that X° = X® = which, in an alternative 
view implies that angular coordinates at T do not shift 
in time as g® A — 0. Our discussion will center around 
the allowed transformations between the induced coordi- 
nates at the worldtube (u, r, and the Bondi system 
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Coordinate and tetrad transformations I 



and 



Let us consider the main gauge freedom admitted in 
our calculation which is of the form 



r 

w(u,(,() 
'1 



o 



r 

,,o\ 



c-c + o 



(30) 
(31) 
(32) 



with a > 0. The possible further transformation of the 
coordinates of the sphere (£, () into itself is not needed 
at this point. 

The condition g ur = 1 in a Bondi system imposes the 
relation 



>J ur 



(33) 



This asymptotic coordinate transformation is associ- 
ated to a corresponding null tetrad transformation; which 
to leading orders is given by 



I = du = adu + a^dC, + a^dC, + O 



O 



(34) 



l A 



(40) 



where 5V is the edth operator of the metric (|28|) . Taking 
into account higher order transformations would include 
an equation of the form 

±=( L A + ho A ); (41) 



l A 



where in principle h could be of order O ( r °). 

The regular dyad at future null infinity in terms of the 
spacetime one, can be given by 



(42) 
(43) 

Then, the transformed regular dyad at future null in- 
finity is given by 



l A = i A . 



5 A = n- 1 5 A 



L 



l A =l A = 



(44) 



(45) 



where we are using f2 



d /i\ i d 
— + o - =-^- 



1 A 
- n + O - 

a \r 



V2P d 
4 

r dC, 



■0 4 



d 



d 



a du d( 
V2P Vwa c Vw 



O 



V 



O 



(35) 



(36) 



since the metric expressed in terms of the new null tetrad 
must coincide with the metric expressed in terms of the 
original null tetrad, it is deduced that 



V 



V Vg° ur 



therefore 



r a 



(37) 



(38) 



The null tetrad transformation equations can be used 
to write the leading order transformation relations for the 
spinor dyad associated to the null tetrad [13, [H|; namely 



4 dya: a 
'a\o A —i A 



(39) 



2. Transformation of ^% 

We can now easily calculate the component \l/4 of the 
Weyl tensor, in leading orders, with respect to the new 
null tetrad, obtaining 

] (46) 



O-iflr ~A~B~C~D 



;,2 4 



(7° 

Hur lT ,0 

a 6 



(47) 



3. Transformations to a Bondi system 

Having analyzed how to transform among different 
frames we can now consider our main task, to relate ^4 
calculated in an arbitrary frame to that which would be 
obtained in a Bondi frame. This frame satisfies, g ur = 1 
and V = 1 which implies 

a = Vg° ur ; a = wg° r (48) 

Which fix both w and a and indicates how different quan- 
tities are to transform. In particular, 

' ' (49) 



^0 = - .jo 

4 v*{gi r y W4 ' 



Thus, knowledge of both g® r and V allows one to obtain 
^4 in terms of the more directly calculated $4. We defer 
to section VI the discussion of how to obtain g^ r and V. 
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Coordinate and tetrad transformations 



metric (the inverse metric), is given by 



After examining the simpler case, we concentrate now 
on the general case, namely we will include the possibility 
of Xq ^ 0. As before, the task at hand is to transform 
to an asymptotic Bondi coordinate (u, f, £, £) and tetrad 
frame (£, h, fh, fh). This transformation is of the form 



1 



u = a(u,(,0 H z h ° ( Z% ) ' 



+ O (r°) 
1 



w(m,C,C) 



C = Co(w,C,C) + o - 



(50) 
(51) 

(52) 



with a > 0. 

If one assumes that £ is an stereographic coordinate of 
the conformal unit sphere (t = const., r = const.) it is 
only necessary to consider an angular transformation of 
the form 



C = Co(«,0 + o(- 

r 



(53) 



The stronger statement is as follows. 

The stereographic coordinates of the sphere can be 
thought of as a map of the extended complex plane. In 
the asymptotic sphere defined by this setting, the angle 
transformation, at u=const., must be conformal. But, a 
map of the extended complex plane onto itself is confor- 
mal if and only if it is a Moebius transformation; that is, 
a transformation of the form 



c 



c( + d' 



(54) 



with ab — be =/= 0. 

We proceed in two stages. 

First, note that the contravariant metric components 
for the standard Bondi like coordinate system is given by 
equations (3.13-18) of [6(; whose inverse is given by equa- 
tions (3.19-24) of the same reference. The only difference 
is that for the general tetrad here considered 
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1 



(55) 



with g ur not necessarily being equal to 1. In order to de- 
duce the required relations we start from the null tetrad 
defined by (TjJJl, 



(D 



.4 



d_ 

dr 



(56) 



9" 
9 V 
9 r 



= 0, 
= A, 

= 0, 
= 2U, 
= X\ 



Then the metric is given by 





1 


Qur 


~ I' 


g rr 


= 0, 


fjri 


= 0, 






Quu 


=-4+ 


9ui 


= -9^ 


9ij 


= {g ij r l 



with i,j,k,l = 2,3, d = det(gij), e,j 
In particular, defining the quantity 



(57) 
(58) 
(59) 
(60) 
(61) 
(62) 



(63) 

(64) 
(65) 

(66) 

(67) 

K"¥); (68) 

-eji and £23 = 1. 



one has that 



1 

w 



(69) 



(70) 



We will be interested in calculating the transformation 
of the radiation field component which depends on 
the transformation properties of the dyad spinor i A . This 
can be read-off from the transformation properties of the 
tetrad vector n. A straightforward calculation reveals 
that the one form n a is 

(n a ) = (g a bn b ) = g ur dr + g uu du + g ui dx l 
+ Ug ru du + X l gi u du + X 1 gijdx J 

= —dr 
A 

-2 j + .V'.V '.,/,., + Uj + .V'.: /,,,.Y'i) du 



+ {-g l3 X j + X j gij) dx l 
= j(dr~ Udu) . 



(71) 



(|2U)l and . Then the contravariant components of the 



In the Bondi frame the following relations are satisfied, 

(l a ) = du, (72) 

(rh a ) = ^e lj CX j du - iey^daj 7 ', (73) 
A A 

(n„) =df- Udu. (74) 
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Additionally, in terms of the original coordinate sys- 
tem, the one-forms associated to the Bondi coordinates, 
to leading order, are given by 



du = a du + ctQ d£ + d£ + O 



(75) 



1 r 

dr = — dr (wdu + w c d( + w^dQ + O(r ), (76) 

w w 2 ^ 



and 



d( = (du + ( c d<; + bd( + 



(77) 



Therefore from the contribution to the metric of £ a fib, 
one deduces that 



a- = q° ■ 

W 



(78) 



that is w is determined by the choice of a and the value 

of g° r - 

In order to study the asymptotic transformation of the 
spinors l A , we examine first the tetrad one- form h a . 



dr — Udu — — dr - (wdu + wrdC + w^d^) + 0(r ) 

w w z ' ^ 



U a du + a.Q d( + d( + O 



--— (dr — ( — w + Ua) du 

r 



(w(d( + w^dC) + (lower orders). 



(79) 



It is important to note that in the general case, the 
vector n a , in the asymptotic region, will have a non 
zero component in the direction of the spacelike direction 
given by the conformal Bondi vectors rh and m. This can 
be seen from the fact that regular tetrad (with hat) at 
future null infinity is related to the standard tetrad in 
the asymptotic region by 



f — f 
m a = Slm Q , 
n a = tt 2 n ai 



(80) 
(81) 
(82) 



one can deduce 



h a (m a ) 



V2P V 



(87) 



where ft = A in this case. It is important here to real- 
ize the following. Extending the hypersurface u =const. 
in the asymptotic region, one defines a section S of fu- 
ture null infinity; in which by construction the vectors 
il~ 1 m a and fi -1 m a are tangent. Also, the regular exten- 
sion of £, namely tt~ 2 £ a , is orthogonal to S. But, since 
future null infinity is a null hypersurface, one has that n a , 
which is a null vector orthogonal to 5*, must be tangent 
to future null infinity. Since that regular extension of the 
Bondi frame to future null infinity is such that fi - rh a 
and fi~ 1 m a will be tangent to S; one deduces that the 
function Co must be chosen such that 



C + A? Cc = 0; 



(88) 



since the vector m a must be orthogonal to n a at future 
null infinity. 

Next, let us study the asymptotic form of n a . Let us 
first note that 



d du d d( d d( d 
du du du du d( du d^ 



(89) 



where the last two terms are different from zero, since 
one has 







d( du d( d( d( 



du du du du d( 
Applying (|8"9")l onto u and using (|90p one can see that 

1 



(90) 



du 

du 



(91) 



and 



dC 

du 



Similarly, one has 



Jp_ _ q,-Jl Co/ 



Co , n 



(92) 



and 



r = rr 2 r, 

m a = n- 1 m a , 
h a = n a - 



(83) 
(84) 
(85) 



where f2 is the conformal factor that defines the regular 
metric g a b = ^ 2 gab- Then, since in the asymptotic region 



(86) 



d(__ 

du 

and thus, 

d_ _ 

du 



Co 



■Coc 



O 



(93) 



'Co c 



-Ss- 



C_d_ _ _C__9 

Cc^c'cc^c 



(94) 
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To leading orders in the asymptotic tetrad transforma- 
tion, n and n are related by 

- d ^f 1 

ou \ r r 



C Co c a <Co< V 

1 / c d ?d 

: -'"-*°V*°* 



Ho< Hoi 



Cc 5 C Cc^' ' 



1 / 1 

= : — n + 0(- 

(95) 

since by virtue of and its complex conjugate, the 
asymptotic terms involving spatial directions cancel. 
Additionally, 

du dudu dC du . du dC ,„„. 
= + -4— = + a c -4; (96) 

9C <9C 9u 9C^C <9C 



thus 



<9£ « du 



(97) 



Last, from 



1 _ (X _ dudC_ + dCd( _ <K _ a_du^ 

~ 8(~ d(du d(d(~ d( C d(~ d( C a c d(' 

(98) 

(99) 



one has 



du 



rt c 



and 



<^C olq C — Q Cc 

ac _ A 
^C C — at 



(100) 



The asymptotic behavior of the vector fh is therefore 

. V2P 8 ( 1 
m = — = + C 



r d{ \r* 
y/2P Vw ( a ( d ad 



1 



O 



w a 

; — r 71 H — ; ~ m 

a Cc ~ a ( C o Cc - C 
1 



° 1 ^ 



(101) 



Since the metric expressed in terms of the new null tetrad 
must coincide with the metric expressed in terms of the 
original null tetrad, it is deduced that 

1 _y_ v\" Cc - c*c Cl _ Vglla Cc - ac CI , (102) 

w a a 2 

where we have used that the tilde system is Bondi like. 
One then can deduce the transformation 



l A 



1 



(103) 



Co 



a — ar~ cka^- 



Consequently the radiation component \&4 transforms 
in the following way 



.f/0 _ o-l,T; ~A~B~C~D 
W 4 — \ l ^ABCDt ILL 



W 



6t-a c 4 !L -ctf-&- 
Q° 

'J it r 



Equivalently, this can be expressed as 



2 4 



2 4 



(104) 



*° 4 = 



Hut ,j,0 

^~2^4 



(105) 



Let us emphasize that Co is chosen to make Xq = 0; 
while a is chosen to make V = 1. 

This complicated transformation of the radiation field 
indicates the convenience of adapting the numerical code 

such that X i = + O(r _1 ); so that C = 0. Otherwise the 
determination of a, needed in (] 1 04[) or (|105[) , would be 
very difficult. 



V. EXAMPLES 

A. Conventions and formulae. Factors of 2 

In section I we presented the standard expressions for 
the radiated four- momentum. Yet, there are differences 
in factors with the definitions employed in numerical 
implementations. The difference arises from the frame 
adopted for the extraction procedure. The analysis pre- 
sented in this work can be used to clarify these differ- 
ences. In a Bondi frame, the expression for the radiated 
four-momentum P a is, 



pa 



1 

47T 



l a a°ir°dS 2 



(106) 



(with ' denoting a derivative with respect to u) This, can 
be re-expressed in terms of ^4 by employing the identity 
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valid in a Bondi frame ^ 



-a . Where we have used 



the "to denote a quantity obtained in a Bondi frame (we 
don't do that for a as here we only use it as obtained in 
the Bondi frame). 



pa 



1 

47T 



I" 



dS 2 



(107) 



In a numerical implementation, as discussed, the read- 
ily available VD^ differs from the Bondi one by 



* 4 



*4 



Thus 



pa = I jo 

47T 



^4 

(a° ) 2 V 3 



du' 



dS 2 



(108) 



(109) 



1. Standard approach, boost ambiguity and resulting factors 

of 2 

We recall that the null tetrads are not uniquely defined 
as they have the boost freedom £ a — > £ a A _1 ,n a — > n a X 
and spin freedom m a — > m a e l7 [2l|. This freedom can 
bring about additional factors in the resulting expression 
for the Weyl scalars, fortunately a Bondi frame naturally 
fixes this ambiguity. In what follows we illustrate how 
this freedom is fixed by ensuring a Bondi tetrad frame 
is adopted in a simple exmaple. This, in passing, will 
make evident the different factors encountered in com- 
monly employed formulae. In the extraction procedure 
it has become custommary to introduce a tetrad frame 
in the following form [42j • First an extraction worldtubc 
is defined as the Cartesian x 2 + y 2 + z 2 = R 2 . Then, 
three vectors (labeled by J) v l j at a hypersurface slice 
are adopted in the following way: 



«i = {-y,x,o),vl = [x,y,z) 1 



det^/V^e^nviv?; 



where 7^ is the induced metric on the spacelike hyper- 
surface at a given time. In our present case, this reduces 
to 



d. 



dl,vl = dl 



(110) 



The next step involves a Gram-Schmidt procedure to con- 
struct three orthonormal vectors v'j with respect to 7^. 
Finally four spacetime vectors are easily constructed by 
r a = '(0,v 2 );0 a = (0,v 3 );<j) a = (0,vi), which together 
with N a (the unit timelike vector normal to the hyper- 
surface) can be employed to construct the tetrad as, 



1 

71 
1 

71 
1 



(N a +r a ) 
(N a - r a ) 



with r 2 = x 2 + y 2 + z 2 . Consider the simplest case of 
a spacetime whose metric, is given by g a b = Vab + h a b 
with h a b a sufficiently fast decaying functions at far dis- 
tances and rjab a flat metric in Cartesian coordinates. In 
this case, to leading order, the tetrad resulting from the 
standard procedure is, 



1 



■Sit) 



d a R ) 



1 



dS- 



(111) 
(112) 
(113) 



rV2 \" a ' sin(0)' 

The main difference with the previous tetrads is that in 
the former one has used a null coordinate u that in this 
case would be of the form u = t — R; with the first null 
tetrad vector given by t a — du a . Therefore, one has 
£' a = -^£ a in this case and so i' A = 2^i A . Then, in 
the calculation of the radiation component of the Weyl 
tensor one would have ^4 = 2^4. The expression for the 
radiated momentum in these coordinates results, 



pa 



1 

167T 



I" 



i$>' 4 du 



dS 2 



(114) 



This then explains the factor of 4 difference in the ex- 
pressions cited in [25|, |28j, |29( with that in eqn (JT^J) (and 
those in the standard literature [l|| |24j, etc.). 



B. Teukolsky waves and wave extraction 

As a second example, we adopt linearized waves on flat 
spacetime as given by the so called Teukolsky waves [43j. 
We adopt this example, coupled with a possible coor- 
dinate transformation to illustrate the differences that 
may arise when suitable contact with a Bondi frame is 
missing. We adopt the simplest expression for the line 
element given by 

(1 + Af rr )dr 2 - 2Bf rg drdO 
(1 + ZCfh, - A)r 2 d0 2 

3C/i -^)r 2 sin(0) 2 # 2 (115) 



ds 2 



dt 2 



1 



where 



and 



frr — 2 

fee - 



A 



3sin(6») 2 ; 
_t\ — f 2 

J d>d> J <hd 



-3 sin(0) cos(0) 
3sin(0) 2 . 



B = - 



C 



4 V 



i?(2) 


3 f(i) 


3F\ 








'^(3) 


3j p(2) 












2F( 3 ) 


9F (2) 









6F 



21F\ 



(116) 
(117) 

(118) 
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with F 



The Riemann tensor for such a line element can be 
straightforwardly constructed. Since to leading order the 
line element is just the flat metric, it is straightforward 
to identify the Bondi frame. A simple calculation gives, 



*° = -sm(6) 2 F^(t~r) 



(119) 



On the other hand, we can calculate '5 4 with the standard 
procedure. We consider two cases: (I) calculating ^ 4 in 
the coordinates (t, r, 8, <fi) and (II) considering the trans- 
formation f — > rg(t) which induces a non-trivial V(t). 
The induced tetrad, to the order that enters in the cal- 
culation of the leading order of ^4 are 

t = ±=(d? + \(l-rg )d% ), 



1 



rgV2 



Taking g = 1 gives the tetrad for case (I). A straightfor- 
ward calculation gives this case, 



#o = £jK«)(t_ r ) sin(6») 2 . 



while for case (II) 



3 F^ (t - rg) sin(6>) 2 

4 9 



(120) 



(121) 



Notice a difference of a factor g 1 appearing in case 
(II) in addition to a factor of 2 when compared to eqn 
(|119|) . These factors result from V — g~ x and g® r = V2g. 
The corrections described in section IV can be used to 
reconcile these differences and obtain Bondi's expression 
in a way mostly independent of gauge issues. In the next 
section we discuss how the different factors involved can 
be obtained in generic settings. 



VI. OBTAINING THE MISSING LINKS 

From our previous discussion, it is clear that one must 
take into account the effects caused by the conformal fac- 
tor V and the induced metric components g ur , g U {. These 
later factors can be read-off from the expression of the 
null one-form h a in the induced null coordinates at T. 
Namely l a = du a defines the function u and the transfor- 
mation from the (t, r, CjC) to the (u, r, £, coordinates 
can be employed to express n a in this system. Then, 
S" 1 =l a dr a ; X* = n a dx\. 

The extraction of the conformal factor V involves more 
work but can be obtained by considering a conformal 
transformation of the sphere and evaluating its curvature 



K = 2 (V 2 + D A D A bxV) 



scalar [4^. Recall that any metric on S 2 is conformally 
related to that of the unit sphere, thus it must be the 
case that g^ = gijR~ 2 = uj~ 2 qij. Thus, the unit sphere 
metric and the induced metric at the worldtube at a con- 
stant time are related by a conformal factor u>. With our 
conventions then V — u. A way to solve for this factor 
is to compute the scalar curvature of the cut S on the 
worldtube. This gives rise to the following relation 

(122) 

where 1Z is the Ricci scalar curvature and Di is the co- 
variant derivative of the metric g^ . Let us observe that 
in two dimensions one has 1Z = 2Ky. 

This equation can be solved to obtain V . Notice that 
this equation admits more than one solution due to the 
rotational group of transformations. One way to solve 
it, and implicitly fix this freedom is motivated in tech- 
niques introduced to find apparent horizons in numerical 
simulation [45l. |46||. In this approach we express V (or a 
related quantity) in terms of a spherical-harmonic expan- 
sion, and then obtain a recursive relation to solve for the 
expansion coefficients. We begin by considering [5l| 

W = V 2 = X lm a lm Y lm (123) 

with Yi m spherical harmonics satisfying V Yi m = + 
^)Yim with V 2 the Laplacian operator with respect to 
the unit sphere metric expressed in the same coordinates 
as those where gij is known (we will refer to this as cjij). 
The solution we seek satisfies H(W) = with 



H = 2W + D l D t InW-TZ. 
Next, we consider the equation 

pH(W) + V 2 W = V 2 W 



(124) 



(125) 

with p a function to be determined in a suitable man- 
ner. The solution we seek turns the equation above into 
an identity. Consider now integrating this equation over 
the sphere having replaced our anzats for W. This will 
provide a recursion relation for obtaining the parameters 
aim (for I ^ 0) as, 



Y lm (pH + V 2 W)dil = -1(1 + l)o, r , 



(126) 



Thus, starting with trial values for {a; m } a new set can 
be obtained through the equation above. Two ingredients 
remain to be provided however, one is how aoo is to be set 
and the other a plausibility argument for the convergence 
of the method. We describe first the latter issue and then 
discuss the former. This argument relies on concentrating 
on the principal part of the system (pH + V 2 W) and to 
exploit p to our advantage. Inspection of these terms 
allows one to re-express them as M lJ DiDjW + S with 



w- 

S = p(^2W 
1 



1 

Tf 



—gV Dl WD 3 W - K 



(127) 



(128) 
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with D the covariant derivative associated with and 
the tensor relating the connections with Di. No- 
tice now that equation (|125j) can be formally solved by 
the iteration scheme 

w (n+i) = (v 2 )" 1 ^ + V 2 )IF (n) (129) 

Here the standard argument applies, the Laplacian op- 
erator smoothes out high frequencies and by properly 
adjusting M y such that the highest order derivatives are 
removed the successive solutions will be smoother which 
is a requirement for convergence. The factor p can be 
chosen so as to cancel particular terms in the second- 
order derivatives, or try to cancel as much as possible 
M lJ (see the discussion in (45|. For instance, if (f- 7 ~ <fi 
the trivial choice p = —1 is ideally suited for this purpose. 

Last we comment on how aoo is obtained. At the first 
step of the iteration, a useful choice is induced by assum- 
ing = al Y 00 . Then, one simply has 

= J (2W (1) - K)Y 0Q <m = 2a} )Q - J TZY 00 dfl (130) 

At subsequent steps, having fixed a™^ 1 for I > 0, equa- 
tion () 1 26(1 can be employed to fix o^q . However, due 
to the non-linearities of the equation this might turn out 
to be difficult. A simpler, more direct route, is to do so 
from 

J [{2W {n+l) ) + D'Di In W {n) - ft) Ybotffl = (131) 

Which can be used to determine a 1 ^ 1 after all other 
coefficients have been obtained. 



VII. FINAL COMMENTS 

In this work we have re-examined the issue of com- 
puting gravitational radiation effects through the use 
of Weyl scalars. The analysis reveals which correcting 
factors are to be accounted for if the coordinates 
adopted do not conform to a Bondi system. To date 
while these corrections have not been taken into account 
in numerical efforts the waveforms obtained both appear 
quite reasonable and, most importantly, agreeing across 
different implementations. What is then the expected 
contribution of the corrections indicated here or the 
relevance of the present discussion in light of these 
observations? 

First, at least the formalism presented in this work 
allows for easily enforcing consistency among different 



implementations. Namely, while coordinates typically 
used vary among different efforts, considering the trans- 
formation to a Bondi system provides a common frame 
for the computation. As a result, comparison of obtained 
signals from different codes would be expedited. Second, 
the correcting terms will have generically non-trivial 
angular dependence. Therefore, the decomposition on 
different multipolar moments would be affected. In the 
case of non-spinning equal mass black holes symmetry 
considerations indicate this might not be a significant 
issue. On the other hand, for different masses and/or 
spinning compact objects the contribution should be 
non-trivial, especially in light of significant kicks ob- 
served which indicate a strong direction dependent of 
the radiated waves. 

It is important to stress here that the corrections are 
case-by-case dependent as they are both sensitive to 
the gauge and initial data employed. Consequently it is 
difficult to assess the role the corrections indicated here 
might have. Nevertheless, consistency with the calcu- 
lation formalism and a simplified frame of comparison 
among different implementations are already key reasons 
for the relevance of the discussion. 

As a side comment we want to stress we have employed 
a convention based on a (—2) signature following [2lLl22l| 
which is the standard signature employed in studies of 
asymptotically flat spacetimes. 

Last, short of considering the corrections discussed 
in this work, as discussed in (47[, one could estimate 
whether these effects may play a role by evaluating the 
norms ||ft — 2||, \\g ur — k\\ (with k = 1 or \/2 depending 
on the boost freedom adopted for l a ) and | \g u A 1 1 - If these 
norms are at the order of the truncation error in a sim- 
ulation, then the correcting factors would certainly not 
be essential. We will examine these issues for different 
binary systems in a forthcoming work[48|. 
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